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Initial magnetic field configurations for 3-dimensional 
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Abstract. We solve, and provide analytical expressions, for current-free magnetic configurations in the context of 
initial setups of 3-dimensional simulations of astrophysical jets involving an accretion disk corona in hydrostatic 
balance around a central object. These configurations which thread through the accretion disk and its corona 
preserve the initial hydrostatic state. This work sets stage for future 3-dimensional jet simulations (including disk 
rotation and mass-load) where launching, acceleration and colhmation mechanisms can be investigated. 
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1. Introduction 

The most successful model of astrophysical jets is that 
of magneto-hydrodynamic (MHD) driven winds as pre- 
sented by Blandford & Payne (1982). In this model, the 
magnetic field is predicted to launch, accelerate and effi- 
ciently collimate the jet. The advent of MHD codes over 
the last decade allowed us to study the details of this 
model through the use of time-dependent MHD simula- 
tions and confirm the importance of the magnetic field. 
One point that has been emphasized in many of these 
simulations is the importance of properly setting up the 
initial state. A numerically stable initial setup makes the 
simulation tractable and allows contact with theory. Here, 
we are particularly interested in simulations as presented 
in Ouyed & Pudritz (1997a, OPI) and Ouyed & Pudritz 
(1997b, OPII) focusing on the technical aspect of the ini- 
tial state. 

The initial conditions as defined in OPI and OPII cor- 
respond to a central object (a proto-star) surrounded by 
a Keplerian disk and an overlying corona in hydrostatic 
equilibrium. The disk has fixed properties and provides the 
boundary conditions for the outflow. The time evolution 
of the initial state is partly depicted by the momentum 
equation (using the standard nomenclature). 



- + (v.V)v 



Vp + pV$ - J X B = 0, 



(1) 



and would, in general, be perturbed by the Lorentz force 
of a magnetic field that threads through the accretion disk 
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and its corona. Therefore, all previous 2-dimensional (2-D) 
simulations (OPI&II) employ initialy current-free, J = 0, 
magnetic fields where the condition for stability becomes 







Vp 



- V$, 



(2) 



at all time. This initial setup was successfully extended 
to 3-dimensional (3-D) simulations in Ouyed, Pudritz, & 
Clarke (2001, OPC). However, the work presented in OPC 
was restricted to the simple case of a uniform magnetic 
configuration where the accretion disks is threaded by ver- 
tical field lines. In this paper, we solve for non-uniform 
current-free magnetic fields in 3-D which would complete 
our 3-D initial hydrostatic corona as defined in (eq (2)). 
The paper is organized as follows: In §2, we describe the 
force-free equations defining our problem. In §3 we solve 
for the current-free case and test the stability of these so- 
lutions by evolving the initial setup in time. We introduce 
and describe the JETSET tool developed for this purpose. 
We conclude in §4. 

2. Force-free configurations 

Force-free configurations are characterized by 

(J X B = 0, V • B = 0), (3) 

corresponding to the standard two cases, J\\B, and J = 0. 
We start with the case where currents are parallel to the 
magnetic field lines. Or, 



B = V X [ — 

a 



(4) 
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(here, a is constant), i.e. the vector potential is given by 
—B (from B = V X A). Thus we first solve for the mag- 
netic field and then calculate the vector potential; recall 
that in deriving B using the vector potential, V • B = 
is numerically guaranteed to within machine round-off er- 
rors. Here and in the rest of the paper, we adopt cylindri- 
cal coordinates (r,0,z) and we assume that the fields are 
separable. 

Using relation (^) and the fact that the magnetic field 
is divergence free we arrive at: 



This is a vector relation and therefore also holds for every 
component. Let us write a component of the magnetic 
field, B, a.s B and making use of our usual assumption of 
separable variables we find: 

B = i?(r)$(0)Z(z), 

R = D{p,\)Jp(V'c 



$ = E{p,X)sm{p4)) 



+ A2rj 
C(p, A) cos(p0). 



(6) 



where A and p are constants of separation, and the un- 
known functions A, C, D, E are to be determined from the 
boundary conditions. The complete solution for B is the 
integral over A ranging from to infinity, and a sum over 
every p (Arfken & Weber 1995): 



B(r,0,z) 



E 

p=0 



D{p, \)Jp + A{p, A) 



X e^^l^l {E{p, A) sin(p<?!)) + C{p, A) cos(p(/i)) 



For every component one finds two constants given 
through Fourier expansion. That is, extra boundary condi- 
tions are required. In our case for example, all components 
of the magnetic field in the accretion disk must be spec- 
ified. However, the integral equations remain difficult to 
treat since the Bessel function is not linearly dependent 
on A, and analytical solutions are not straightforward. We 
thus turn to the current-free case (a = 0) where the equa- 
tions can be simplified. 

3. Current-free configurations 

Here, J = is a necessary and sufficient condition to 
guaranty the existence of a scalar field (ip) with the prop- 
erty, Vv? = B. As can be seen in appendix ^ this greatly 
simplifies our task of finding the initial configurations. We 
investigate two kind of current-free magnetic fields, with 
and without toroidal component (i?^). 

3.1. B^^O 

Here we do not demand that the toroidal field is zero in 
the corona. We use the following boundary condition: 



</>,o 



{b r^^i) X sin(fc0). 



(7) 



where B^.o is the toroidal magnetic field in the disk, and b 
a normalization factor. The general solution for the scalar 
field is given in appendix ^ where we also explain the 
choice of such a (j) dependence. We find, 

ip{r,4>,z) — q{p,,k)bcos{k4>) 



+ \z\ 



^,'k-p + lu + k, 



(5) q{p,k) = 2^-^-+ 



^^^r{}{k + p + 2))r{k-p + 2) 



r(i(fc-^))r(fc + i) ' 

for 1 - < < i , (8) 

where 2^1 is a hyper geometric function and F the gamma 
function. The k = 2 case, for example, is a simple solu- 
tion. The corresponding configurations once implemented 
in the simulations remain stable in time. In the B^ ^ 
case, however, the corresponding jet simulations (in- 
cluding mass-load and disk rotation) are prone to pinch 
forces within few disk rotations. That is, the resulting 
dynamics is more complex, and the jet more difficult to 
track/investigate numerically. The B^ = configurations, 
also adopted in OPI&II, turned out to be useful in many 
ways; it allows for instance to demonstrate that the jet col- 
limation can result from the self-generated toroidal field. 
These we consider in details next. 



3.2. B^ = 

In the current-free initial setup the differential equation 
for the scalar field of the magnetic field {(p) reduces to the 
Laplace equation (since V • B = 0): 



V2^(r,(/),z) = 0. 



(9) 



The corresponding scalar potential is found to be 
(Appendix]^, 



(p{r,<l),z) = 2'^6z^ 



F(l-hf)F(-M) ^A-p p 



F(l-f) 
1 



(10) 



Solution above corresponds to the following boundary con- 
dition: 



H,.n = b 



(11) 



where -6^,0 is the radial component of the magnetic field 
in the accretion disk. Recall that we do not simulate the 
accretion disk (a fixed boundary in all of our simulations) , 
but only the jet and the overlying corona (z > 0). 

We examined configurations for 2 different values of 
p (the radial dependence of the magnetic field); an open 
{p = —1) and a closed {p — —2) configuration. The closed 
configuration has the same radial dependence as a dipole 
field, which is expected (in first approximation) around 
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Fig. 1. Current-free solutions (B^ = 0) - the open case 

(/i = —1): Shown here are the magnetic field lines at r = 
(top) and at T = 100 (bottom). In this figure and the rest 
of figures, the blob in the center represents an iso-surface 
density around the central object (the proto-star is less 
than few pixels in size and cannot be seen in the figure). 
The rotation axis is shown crossing the iso-surface density 
into the central object and is in the plane of the paper. 
The accretion disk not shown here is to the far left and is 
perpendicular to the rotation axis. 



stars however ours is strictly a solution of the disk bound- 
ary condition. The open configuration is the 3-D ana- 
logue/extension of the open configuration used in OPI. 
We implemented these 2 configurations in our initial setup 
and let them evolve in time after applying a small per- 
turbation to the density. No mass injection and no disk 
rotation were set, and = Vr = = in the corona 
{z > 0.0). For the simulations we used the time-explicit 
Eulerian MHD-codc ZEUS3D (Stone & Norman 1992). 

Cartesian coordinates, (x, y, z), are used for all simula- 
tions. While being the natural system to use to avoid any 
directional biases, it docs introduce some of its own prob- 
lems not encountered in the 2-D cylindrically symmetric 
simulations. We refer the interested reader to OPC for the 
technical reasons underlying the choice of Cartesian coor- 
dinates for such simulations. The disc is taken to lie along 
the x-y plane, and the disc axis corresponds to the .z-axis. 
In units of the inside radius of the disc, r^, the simulated 
region has dimensions (—15 : -|-15, — 15 : -t-15,0 : +60), 
and is divided into (95, 95, 120) uniform rectangular zones. 

Fig.l shows the magnetic field lines for the open config- 
uration at time t = (in units of inner Kepler period) and 



Open 




Fig. 2. Current-free solutions {B^ = 0) - the open 
case (^ = — 1; top) and the dipole case (/x = —2; bottom). 
Shown evolving in time are the maximum values for the 
velocity (solid line), the magnetic field strength (dotted 
line), and the density (dashed lines). The dot-dashed line 
shows the maximum speed for the propagation of Alfven 
waves. 
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Fig. 3. Current-free solutions (i?^ = 0) - the dipole 
case (/i = — 2). Shown are the magnetic field lines at r = 
(top) and at T = 100 (bottom). 



at time r = 100. The magnetic configuration changes only 
slightly in time and the generated currents are negligible. 
This is further demonstrated in Fig. 2 (top panel) where 
we show the time evolution of the maximum values for 
velocity, density and magnetic field strength. The induced 
velocity and slight variations in the magnetic field topol- 
ogy arc signature of Alfven waves which as expected have 
no dynamical effects on the coronal material. The solution 
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stays force-free up to r = 300 and no motion/jet is gen- 
erated - as can further be checked from the movies of the 
simulations (at http://www.fys.ku.dk/"" svart/Jets/). The 



closed configuration (lower panel in Fig. 2) is also stable in 
time although the induced errors are 20% larger than in 
the open case (Fig. 3). Here as well the initial hydrostatic 
balance is preserved up to r = 300. 

The simulations described above are easily imple- 
mented numerically and if it is left unperturbed, the 
corona will remain in perfect numerical balance to within 
machine round-off errors. In any case, when disk rotation 
and mass- load are taken into account, the small pertur- 
bations are quickly and completely washed out by the jet 
dynamics (r 50). The complete set of simulations with 
and without disk rota tion and mass-loading can be visual- 
ized and compared at http://www.fys.ku. dk/^ svart/Jets/ . 



3.3. The JETSET tool 

We developed a tool, namely JETSET, that generates ini- 
tial states as described above. The JETSET main frame 
is shown in Fig. 4. Once the grid dimensions, the phys- 
ical scales and the appropriate coronal and disk param- 
eters have been specified, JETSET performs a Newton- 
Raphson method to find the correct density distribution 
while the corresponding magnetic field configuration is 
computed using the approach described in previous sec- 
tions. The resulting Data (density, specific energies, ve- 
locity, and magnetic field) describing the initial setup is 
stored in an HDF (Hierarchical Data Format) file which 
can then be read by the user's code (such as ZEUS). 
The magnetic field lines can be visualized (see Fig. 4) 
as well as the density distribution around the central ob- 
ject (see Fig. 5). JETSET is available (down-loadable) at 



http://www.nordita.dk/"" ouyed/ JETTOOLS/ . Included in 



the package are README and HELP files. 
4. Conclusion 

In this paper, we solved analytically for force- free solutions 
of magnetic configurations which can be implemented in 
3-D simulations of astrophysical jets (disk winds). These 
configurations which thread the accretion disk and the 
corona, we showed, do not perturb the initial hydrostatic 
balance and are stable in time. While idealistic (developed 
for ease of implementation and computation), they con- 
stitute the first stage towards testing the effects of differ- 
ent magnetic configurations on the simulated jets in 3-D. 
Realistic configurations ought to reproduce basic features 
of astrophysical jets, such as their cylindrical shape, their 
knotty structure, and their stability. 

Acknowledgements. We thank R. E. Pudritz and C. Rogers for 
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Appendix A: Current-free configurations 

A.l. B^^O 

Here, the problem reduces to finding a general solution to 
Laplace's equation: 



VV = 0. 



(A.l) 



-Al 



Making a separation of variables in cylindrical coordi- 
nates, 

^ = i?^(r)$^(0)Z^(z) , (A.2) 
we find (Gradshteyn & Ryznik), 

^(r, 0,z) = / V D{p, X)Jp{Xr)A{p, X)e 

p=0 

{E{p, X) sm{p(f>) + C{p, X) cos(p0)) dA, (A.3) 

where A and p are constants of separation, and Jp(Xr) 
the Bessel function. The unknown functions A,C,D,E 
which can in fact be combined in two constants (see be- 
low) within the Fourier space are to be determined from 
the boundary conditions. Equation ( A.3 ) is a superposi- 



tion of the solutions to R^p, and Z^,. It shows that the 
constants may be determined from a single boundary con- 
dition, for example the configuration of the toroidal mag- 
netic field in the disk B^ q = B^{r, (f>, 0). 

We now look at the i?^ = 'Vip\^ component: 



cos(p0) ( p I 5i,p(A)e-^l^l Jp(Ar)dA ) + 



sm{pcl>) (-P / 52,p(A)e-^l^ljp(Ar)dA 



p=0 



where Si p{X) = D{p,X)E{p,X)A{p,X) and S2 p{X) = 
D{pA)C{p,X)A{p,X) 

Equation above is simply a Fourier expansion in the 
toroidal dependence of the magnetic field. The integrals 
before the trigonometric functions are the Fourier coeffi- 
cients; so we may write : 

/ 5i.p(A)Ae-^l^ljp(Ar)dA = — / sin(p,/))d(^, 
Jo T^P Jo 

52,p(A)Ae-^ 



I r'' 

Jp(Xr)dX = — / B^cos{p(j))d(j), 
T^P Jo 

where the right hand side is the expressions for the Fourier 
coefficients. Note that S'i.p(A) and 5*2, p(A) are the Hankel 
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Fig. 4. JETSET TOOL: JETSET sets initial states for astrophysical jet simulations in 3-D. The initial set up 
consists of a corona in hydrostatic balance around a central object with current-free magnetic fields threading through 
the corona and the underlying accretion disk (fixed boundary). JETSET solves for the correct density distribution 
around the central object and the appropriate current-free magnetic field configuration once the two parameters, k 
and fi, are specified (see text). The coronal material might consist of one pressure component (71 = 5/3, 72 = 0) 
or two-pressure component (71 7^ 0, 72 7^ 0, where 71 and 72 are the corresponding adiabatic indices). The initial 
magnetic field configuration can be visualized and saved into a file. Figures similar to the upper panels in Fig. 1 and 
Fig. 2 can be generated (and saved into files) by JETSET as evident from the renderer to the right. Further details on 
JETSET can be found in the README and HELP files included in the JETSET package which can be down-loaded 
at 



http://www.nordita.dk/'^ ouyed/JETTOOLS/ . 




Fig. 5. JETSET TOOL - Density contours as generated by JETSET for a single pressure component corona with 
71 = 5/3. The left panel is in a plane parallel to the disk surface while the panel to the right shows the density contours 
in a plane perpendicular to the disk surface and containing the disk rotational axis. 
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transforms of the Fourier components of B^{r, cj), 0). That 
is, for any given toroidal configuration of the magnetic 
field in the disk, Si^p{X) and 52,p(A) are given by equations 
above thus providing a simple procedure for determining 
the scalar field. 

Finally, note that the simplest magnetic configuration 
possible is given by Sp{X) = OVp ^ fc which gives us 
the boundary condition investigated in this paper -B^.o = 
(b r'^^-^) X sin(fc(/)) (including the fc = case with no 
toroidal dependence); 6 is a scaling factor. 



A.2. B^ = 

The corresponding scalar field is obtained by cancelling 



the (j) dependence in (A.3). The components Br and 
are then direct derivative of ip. For example, one can show 
that 

B,= 5(A)AJo(Ar)e-^l"ldA , 

"'0 



where S'(A) is the Hankel transform of i?z(r, (f), 0). 



